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ABSTRACT 

We investigate whether stehar-mass black holes have to receive natal kicks in order 
to explain the observed distribution of low-mass X-ray binaries containing black holes 
within our Galaxy. Such binaries are the product of binary evolution, where the mas- 
sive primary has exploded forming a stellar-mass black hole, probably after a common 
envelope phase where the system contracted down to separations of order 10 — 30 R©. 
We perform population synthesis calculations of these binaries, applying both kicks 
due to supernova mass-loss and natal kicks to the newly-formed black hole. We then 
integrate the trajectories of the binary systems within the Galactic potential. We find 
that natal kicks are in fact necessary to reach the large distances above the Galactic 
plane achieved by some binaries. Further, we find that the distribution of natal kicks 
would seem to be similar to that of neutron stars, rather than one where the kick 
velocities are reduced by the ratio of black hole to neutron-star mass (i.e. where the 
kicks have the same momentum). This result is somewhat surprising; in many pictures 
of stellar-mass black-hole formation, one might have expected black holes to receive 
kicks having the same momentum (rather than the same speed) as those given to 
neutron stars. 

Key words: X-rays: binaries - stars: neutron - supernovae: general - Galaxy: dy- 
namics - binaries: general - black hole physics 



1 INTRODUCTION 

It has long been known that neutron stars receive kicks at 
birth in the range ~ 200 — 400 km/s (so called natal kicks), 
when they are formed in core-collapse su pernovae, for exam- 
ple via proper motion studies of pulsars (jCordes et al.lll993l . 
iLvne fc Lorimeii [l99J). Whether stellar- mass black holes 
(for brevity, black holes hereafter) receive these kicks too 
is still a matter of debate. Black holes can be studied via in- 
teracting X-ray binaries which contain them. There are sev- 
eral known X-ray binaries which are kn own to contain black 
holes or contain black hole candidates (jJonker fc Nelemand 
I2OO4I . lOzel et al.ll2010l ). In these systems, the massive pri- 
mary has evolved to form a black hole via a core-collapse 
supernova and material is currently flowing from the lower- 
mass secondary (typically via Roche-lobe overflow) onto the 
black hole via an accretion disc (f or a detailed review on the 
evolu tion of compact binaries see iTauris fc van den Heuvej 
120061 ) . When the primary explodes as a supernova, the mass 
loss from the system can unbind the binary or at least give 
it a kick (as the mass lost has a net momentum in the rest 
frame of the binary). In addition, any natal kick received by 
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the black hole will affect the orbital properties of the binary 
and its orbit within the Galaxy. By studying the orbit of 
a binary, or even its location within the Galaxy, one might 
obtain a limit on the range of allowed natal kicks. A number 
of studies have considered the motion of individual binaries 
within the Galaxy. 

iBrandt et all Il995l considered GRO J 1655-40 (Nova 
Sco) and concluded that a natal kick more easily accounted 
for the high space velocity of the binary. 

iNelemans et al.lll999l studied Cygnus X-1 and conclud- 
ing that a natal kick was not necessary to explain its space 
velocity. 

IWillems et ahlbOOSl considered GRO J1655-40 and sug- 
gested that although a natal kick is not formally required 
to produce the system as observed today, the inclusion of a 
(modest) natal kick more readily explains the system. They 
also placed an upper limit on the natal kick of ~ 210 km/s. 

iDhawan et"ai]|2007l considered GRS 1915+105. They 
concluded that any peculiar motion of the binary was more 
likely due to later scattering within the Galactic disc than 
a natal kick when the black hole formed. 

The binary XTE JU 18-1-480 is located at a very 
high latitude ( 1.5 k pc above the Galactic disc, see 
iRemillard et al.l [2000|) and it has a high space velocity 
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l|Mirabel et al■ll2001^ . [Gualandris et al.ll2005l concluded that 
for this system a black ho le natal kick was required. More 
recently. iFragos et al.ll2009l placed the value of the natal kick 
in t he rauRe 80 - 310 km/s. 

IWong et alj I2OI0I considered Cygnus X-1. They found 
that in this case the black hole progenitor could have re- 
ceived a relatively small natal kick (few tens of km/s with 
an upper limit of 77 km/s). If the sy stem originated in the 
Cyg OB3 association (jMirabel fc Rodrigues 2003 ') , then the 
upper limit on a kick is reduced to 24 km/s. 

In this Paper, we consider the population of black 
ho le X-ray binaries as a who l e (following the approach 
of IWhite fc van Paradi"i3 1 19961 . iJonker k. Nelemanj |2004 
IZuo et al.ll2008l l rather than consider the kinematics of an 
individual system. We synthesize a population of black-hole 
low-mass X-ray binaries (BH-LMXBs) , using various natal 
kick distributions, and integrating the systems within the 
Galactic potential. We then compare their locations within 
the Galaxy to a catalogue of kn own black-hole X -ray bina- 
ries having measured distances l|Ozel et al.ll20ld ). 

The paper is arranged as follows. We review the current 
state of observations of X-ray binaries containing either neu- 
tron stars or black holes in Section (2] Our treatment of the 
motion of stars in the Galactic potential is given in Section 
O In Section[3]we review the effects of both natal and super- 
novae mass-loss kicks on binaries. In Section[5]we present the 
results of the binary population synthesis which we discuss 
in Section [S] The paper is concluded in Section [T] 
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Figure 1. Galactic distribution of NS-LMXBs (open circles) and 
BH-LMXBs (filled circles). The radial distance from the Galactic 
centre is i? = ^Jx^^+y^ , where x and y are the cartesian coordi- 
nates in the Galactic plane, whereas the distance from the plane 
of the Galaxy is 2 = dsinb. Two neutron star binaries fall off the 
figure: XTEJ2123-058 and CygX-2. Solid lines account for the un- 
certainty of the distance from the Sun for each binary. The Sun 
is indicated as a square. 



2 THE OBSERVED BINARIES 

In our Galaxy there are 16 dynamically confirmed black 
holes in LMXBs and 33 NS-LMXBs, wh ose distance and 
Galactic position is known ; see respectivelv lOzel et al.ll2010l . 
Ijonker fc Nelemansll2004l . and references therein (in partic- 
ular, Jonker & Nelemans consider only NSs not found in 
globular clusters). We present the binaries in Tables [1] and 
[51 along with their angular distribution, their distance from 
the Sun and their position, both in Galactic coordinates and 
in cylindrical ones {R refers to the radial distance from the 
Galactic centre and z refers to the distance from the Galac- 
tic plane). Concerning the black hole binaries, uncertainty 
on the distance is taken from Ozel et al.; for the neutron star 
binaries, Jonker & Nelemans calculated the distance assum- 
ing two different Eddington peak fluxes, getting a maximum 
and a minimum value for the distance, of which we take the 
median value. 

Using the values in Table [l] and Table [2j we plot the 
Galactic distribution of the binaries in Figure [l] In repre- 
senting the observed system on the [R, z) plane we propa- 
gate the uncertainty on the distance into an uncertainty on 
R and on z, the corresponding range of values for R and z 
being represented as a solid line. The z-distribution of BH- 
LMXBs appea rs similar to the NS-LMX Bs one (as already 
pointed out bv lJonker fc Nelemang|2004l . who calculated the 
rms z for the two samples); the Kolmogorov-Smirnov prob- 
ability P for the two distributions to be the same is indeed 
convincing: P ~ 0.81 (that increases to 0.85 if binaries lo- 
cated at 2: > 2.0 kpc are excluded from the test). Concerning 
the R distribution, we may observe that NS-LMXBs seem to 
be more concentrated towards the Galactic centre with re- 



spect to the BH-LMXBs. Th is is very likely an observa tional 
bias, as already suggested bv ljonker fc Neleman32004l . since 
the BH binaries that we consider are those for which a dy- 
namical measurement of the BH mass exists, and this biases 
the binaries towards those closer to us. 



3 INTEGRATING ORBITS WITHIN THE 
GALAXY 

In this section, we will give a short overview of how the or- 
bital trajectories within the Galaxy are calculated. For the 
Galactic potent ial, we make use of the model proposed by 
IPaczvnskill 199d . The potential is modeled as the superposi- 
tion of three components due to the disc ($d), the spheroid 
(<&s) and the halo ("^h), as given below 

$d(i?, z) = ^^^ = (1) 



/i?2 + (ad + y/z^ + fe|) 
where a^ = 3.7 kpc, &d = 0.20 kpc, and Md = 8.07x lO^^Mg. 



$s(i?,z) = 



GM, 



where as 
1O^°M0. 



' R^ + [os + V^n~bT} 
0.0 kpc, fes = 0.277 kpc, and Ms 



$h(r-) = 



GAfc 



1 , / , r~ \ Tc I r 

— m 1 H — T H arctan — 



(2) 
1.12 X 

(3) 
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Table 1. Observed properties of BH-LMXBs. 



Name 



(dog) 



b 

(deg) 



d 
(kpc) 



Ad 
(kpc) 



R 

(kpc) 



(kpc) 



Ref. 

(distance) 



4U1543-47 


330.0 


+5A 


7.5 


0.5 


3.92 


0.70 


[1] 


XTEJ1550-564 


325.9 


-1.8 


4.4 


0.5 


5.0 


-0.14 


[2] 


GROJ1655-40 


345.0 


+2.5 


3.2 


0.5 


4.98 


0.13 


[3] 


1659-487 


338.9 


-4.3 


9.0 


3.0 


3.25 


-0.67 


[4] 


1819.3-2525 


6.8 


-4.8 


9.9 


2.4 


2.14 


-0.82 


[5] 


GRS1915-I-105 


45.4 


-0.2 


9.0 


3.0 


6.62 


-0.03 


[6] 


GS2023-I-338 


73.1 


-2.1 


2.39 


0.14 


7.65 


-0.09 


[7] 


GROJ0422-f32 


166.0 


-12.0 


2.0 


1.0 


9.91 


-0.41 


[8] 


A0620-003 


210.0 


-6.5 


1.06 


0.12 


8.92 


-0.12 


[9] 


GRS1009-45 


275.9 


-1-9.4 


3.82 


0.27 


8.48 


0.62 


[10] 


XTEJ1118-I-480 


157.6 


+62.3 


1.7 


0.1 


8.73 


1.50 


[11] 


1124-683 


295.3 


-7.1 


5.89 


0.26 


7.63 


-0.73 


[10] 


XTEJ1650-500 


336.7 


-3.4 


2.6 


0.7 


5.71 


-0.15 


[12] 


1705-250 


358.2 


+9.1 


8.6 


2.1 


0.55 


1.36 


[13] 


XTEJ1859-I-226 


54.1 


+8.6 


8.0 


3.0 


7.23 


1.20 


[10] 


GS2000-I-251 


63.4 


-3.0 


2.7 


0.7 


7.21 


-0.14 


[13] 



References ffrom | 62el et al.ll2010 M: [1] Orosz 2010 priv ate commu nication, [2lJOrosz_^t_al 



2010, [3] Hi ellming fc Rupcn 1995, [4] Hyncs 2005, f5] |Orosz et"a l. 2001, [ 61 iFender et al 
~999 . [71 .Miller- Jones et al..2009. . [81 .Webb et al 20001 . [91 ICantrell ct al.l20ld . [lOl lHvnei 
20051. [ll] lGelino et al.ll2006l . Tl2] iHoman et al.lll999r [13l lBarret et al.lll996l . 



where re = 6.0 kpc and Mc = 5.0 x IO^^Mq. 

When integrating the trajectory of a binary within the 
Galaxy, we make use of the cylindrical symmetry of the po- 
tential. The equations of motion which are thus integrated 
are given below 



dR 

11 



= Vr, 



dvR 
dt 



dR 



+ ii 

R^ 



dz 
dt 



dv^ 
~dt 



9£ 

dz 



(4) 



(5) 



where R and z are the cylindrical coordinates of the binary, 
jz is the z-component of the angular niomentum of the bi- 
nary, and $ = $d + *a + $h. 

It will turn out that typical kick velocities that the bi- 
nary receives when the primary explodes as a supernova 
are comparable to the circular orbital speed in the Galaxy 
(~ 200 km/s). This implies that kicks can significantly af- 
fect the trajectory of the binary within the Galaxy. We can 
get an idea of the maximum z reached by the binary as 
a function of the initial peculiar velocity. We integrate the 
equations of motion of the binary for ~ 10 Gyrs (which is 
the typical main sequence MS-time of a sun- like star), using 
a 4th-order Runge-Kutta integrator developed by SR, and 
assuming the system was born right in the Galactic plane 
with a peculiar velocity perpendicular to the plane v^. We 
perform the integration for different values of the velocity 
(nx = 20, 40, 100, 200 km/s) and of the initial position Rt=o 
over the plane, writing down the maximum z reached over 
the trajectory (see Figure [2|. We see how Zmax is a rather 
strong function of the initial position: for a fixed value of 
v_i_, 2max gcts Smaller as the binary gets deeper in the po- 
tential well (i.e. smaller values of Rt=o). It is clear from this 
figure that binary kick speeds in excess of 200 km/s will be 







4 6 
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Figure 2. Maximum z reached by a binary over its trajectory. 
The object has been kicked perpendicularly to the Galactic plane, 
for four different magnitudes of the kick (v± = 20, 40, 100, 200 
km/s) . 



required in at least some of the observed systems shown in 
Figure 1. 
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Table 2. Observed properties of NS-LMXBs. 



Name 


1 


b 


d 


Ad R 


z Ref. 




(deg) 


(dog) 


(kpc) 


(kpc) (kpc) 


(kpc) (distance) 


EXO0748-676 


279.98 


-19.81 


7.95 


2.3 9.96 


-2.70 [1] 


2S0918-54 


275.85 


-3.84 


5.05 


1.5 9.0 


-0.34 [2] 


Cir X-1 


322.12 


0.04 


9.15 


2.7 5.67 


0.00 [3] 


4U1608-522 


330.93 


-0.85 


3.3 


1.0 5.36 


-0.04 [4] 


Sco X-1 


350.09 


23.78 


2.8 


0.3 5.49 


1.13 [5] 


4U1636-53 


332.91 


-4.82 


4.3 


1.2 4.62 


-0.36 [6] 


4U1658-298 


353.83 


7.27 


9.85 


2.9 2.01 


1.25 [7] 


4U1702-429 


343.89 


-1.32 


6.2 


1.8 2.67 


-0.14 [8] 


4U 1705-44 


343.32 


-2.34 


8.4 


2.4 2.40 


-0.34 [8] 


XTEJ1710-281 


356.36 


6.92 


17.3 


5.0 9.20 


2.08 [8] 


SAXJ1712. 6-3739 


348.93 


0.93 


6.9 


2.0 1.81 


0.11 [9] 


1H1715-321 


354.13 


3.06 


6.0 


1.8 2.13 


0.32 [10] 


RXJ1718.4-4029 


347.28 


-1.65 


7.5 


2.2 1.79 


-0.21 [11] 


4U1728-34 


354.30 


-0.15 


5.3 


1.6 2.78 


-0.01 [12] 


KS1731-260 


1.07 


3.65 


6.2 


1.8 1.81 


0.39 [13] 


4U 1735-44 


346.05 


-6.99 


9.4 


2.8 2.48 


-1.14 [8] 


GRS1741. 9-2853 


359.96 


0.13 


7.75 


2.3 0.25 


0.01 [14] 


2E1742. 9-2929 


359.56 


-0.39 


8.05 


2.3 0.08 


-0.05 [8] 


SAXJ1747.0-2853 


0.21 


-0.24 


8.75 


2.5 0.75 


-0.04 [15] 


GX3+1 


2.29 


0.79 


5.05 


1.5 2.96 


0.07 [16] 


SAXJ1750.8-2900 


0.45 


-0.95 


6.1 


1.8 1.90 


-0.10 [17] 


SAXJ1752. 3-3138 


358.44 


-2.64 


9.25 


2.7 1.26 


-0.43 [18] 


SAXJ1808. 4-3658 


355.38 


-8.15 


3.15 


0.9 4.90 


-0.45 [19] 


SAXJ1810.8-2609 


5.20 


-3.43 


5.95 


1.7 2.15 


-0.36 [20] 


4U1812-12 


18.06 


2.38 


4.0 


1.2 4.38 


0.17 [21] 


XTEJ1814-338 


358.75 


-7.59 


9.6 


2.8 1.53 


-1.27 [22] 


GX17-I-2 


16.43 


1.28 


13.95 


4.1 6.67 


0.31 [23] 


SerX-1 


36.12 


4.84 


11.1 


3.2 6.59 


0.94 [8] 


AqlX-1 


35.72 


-4.14 


5.15 


1.5 4.86 


-0.37 [24] 


4U1857-I-01 


35.02 


-3.71 


8.75 


2.5 5.08 


-0.57 [25] 


4U1916-053 


31.36 


-8.46 


8.8 


2.6 4.56 


-1.29 [26] 


XTEJ2123-058 


46.48 


-36.20 


18.35 


5.3 10.96 


-10.83 [27] 


Cyg X-2 


87.33 


-11.32 


13.35 


3.9 15.02 


-2.62 [28] 
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4 KICKS RECEIVED BY SURVIVING 
BINARIES 

In this section we consider the effects of the supernova explo- 
sion on the binary. We will see how the rapid mass loss from 
the supernova alone could impart a kick on some systems 
whilst breaking others up. In addition, any kick imparted to 
the neutron star or black hole on its formation (i.e. a natal 
kick) will also play a role, both in adding to the overall kick 
received by the binary, and in some cases ensuring that the 
binary remains bound. 

It is important to note that a conspiracy of three ve- 
locities will have an important role; namely the coincidence 
that the following three speeds are comparable: the speed 
of a circular orbit in the Galaxy, the typical orbital speed 
within a tight stellar binary when the primary explodes as a 



supernova, and the characteristic kick speed the binary re- 
ceives. This coincidence implies that kicks will significantly 
affect the orbit of the binary within the Galaxy. 

We begin by considering the case of zero natal kick. In 
other words, where any kick is due solely to the rapid mass 
loss occurring during the supernova explosion. We will refer 
to this as the meiss-loss kick Knik (also called Blaauw kick, 
lBlaauwlll961f ), which is given by the expression below 



Kx, 



AM Ma GM 



M' M 



(6) 



where M is the total mass of the binary at the point of 
the supernova explosion, M' is the total mass of the binary 
after the supernova explosion, AM is the mass lost during 
the supernova explosion (i.e. AA/ = M — M'), M2 is the 
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mass of the secondary, and a is the binary semi-major axis 
at the moment of the supernova explosion. 

Knik can't be too large, since the binary must remain 
bound after the supernova: the mass loss must be less than 
half of the initial mass. If we agree on a common envelope 
phase having shrunk the binary down to an orbital sepa- 
ration of ~ 10 R0, the resulting typical mass loss kicks for 
BH-LMXBs are of the order of 20-40 km/s (for NS-LMXBs 
they are typically higher). Looking at figure [21 we immedi- 
ately realize how kicks of this size cannot make the highest-z 
black hole binaries: in the optimal case of Knik — 40 km/s 
perpendicular to the Galactic plane, the maximum z reached 
over the trajectory never exceeds 1 kpc (however, we do see 
binaries in the halo of our Galaxy at larger values of z, see 
table [TJ. 

We consider now the case where the neutron star or 
black hole produced in the supernova receives a natal kick. 
If we assume the orientation of the natal kick is random with 
respect to the orbital plane, the natal kick Vlik combines with 
the mass loss kick Knik as given below: 



Vk 



M' 



^ ,,, V^nk.xVmlk 



(7) 



where we have chosen the x axis aligned with the orbital 
speed of the BH progenitor and the y axis along the line 
connecting the two stars at the moment of the supernova 
explosion. 

Many distributions to model neutron star natal kicks 
have been proposed. For example, [Hansen fc Phinnevlll997l 
modeled the natal kick as a Maxwellian distribution peaked 
at 300 km/s. To solve the retention problem in globular 
clusters as well as the low eccentricity of a subclass of Be 
X-ray binaries, two-peak distributions have also been pro- 
posed where one peak occurs a.t a som ewhat lower velocity 
i|Pfahl et al.| [2002a. Pf ahl et al.|[2002bh . 

We consider two different natal kick distributions here: 
one is the Hansen & Phinney di stribution, the other is a 
bimodal distribution proposed bv lArzoumanian et al.ll2002l 
which has a lower peak at ~ 100 km/s and the higher 
peak at ~ 700 km/s. We also consider modified versions 
of the above two distributions, which we term momentum- 
conserving kicks MCKs, where we assume that the momen- 
tum imparted on a black hole is the same as the momentum 
given to a neutron star using the two distributions. Thus the 
kick velocities will be reduced: T4k,bh = (Mns/Mbh)Vnk,ns. 
For example, a 7 M© black hole receives a natal kick reduced 
by a factor of 5: for a neutron-star natal kick of 300 km/s, 
the black hole would receive a smaller kick of only 60 km/s. 
We show in figure [3] the natal kick distributions which we 
use. 

It is important to recall that a large fraction of binaries 
are broken up when the primary explodes as a supernova. 

Considering a population of binaries where M\ = 
11 M0, M2 = 1.5 Mq, Mih = 7.8 M0, a = 10 R0 (Mi 
is the mass of the progenitor of the black hole, M2 is the 
companion star and a is the pre-SN orbital separation), we 
impart the BH of each binary a NK drawn randomly from 
each of our four distributions. The fraction of systems re- 
maining bound for each of the kick distributions is shown 
in Table (3] We also include the case where a 1.4 Mq neu- 
tron star is produced instead of a black hole (taking as the 
progenitor mass 3.5 Mq). A larger fraction of binaries re- 
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V„, (km/s) 

Figure 3. Natal Kick distributions used in our binary population 
synthesis calculations. Solid line corresponds to Arzoumanian dis- 
tribution, dotted line to Hansen & Phinney and the two dashed 
lines to these two distributions but with kick speeds reduced, as- 
suming that the momentum imparted to the black hole is the 
same as the momentum imparted to the neutron star. 



Table 3. Fraction of systems that stay bound after the SN. 



Fraction of bound systems 




BH 


NS 


Hansen 


58% 


30% 




56%" 


35%" 




64%'' 


3%' 


Bimodal 


54% 


29% 




56%" 


30%" 




52%'' 


10%'' 


Hansen MCK 


99% 


- 


Bimodal MCK 


95% 


- 



" For a NK lying in the orbital plane. 

*" For a NK perpendicular to the orbital plane. 



main bound for binaries containing black holes rather than 
neutron stars owing to the greater binding mass. 

We show in figure [4] the distribution of kick velocities 
for BH-LMXBs that we obtain drawing from each of the 
four natal kick distributions. One should in particular note 
how the kick velocities for the momentum conserving kicks 
are typically lower than ~ 100 km/s. 



5 BINARY POPULATION SYNTHESIS 

In this section we discuss the calculation of the synthetic 
population of black hole low-mass X-ray binaries. We pro- 
duce a population of BH-LMXBs considering their forma- 
tion within the Galactic disc at a range of radii. For each 
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Figure 4. Peculiar velocity gained by the binary after an asym- 
metric supernova. Natal kicks have been drawn from Arzouma- 
nian distribution (solid line), Hansen & Phinney distribution 
(dotted line), whereas the dashed lines correspond to the MCKs. 
The total number of systems for each curve has been normal- 
ized to 1000 and only systems that stay bound after the SN are 
represented. 



binary, we randomly draw the black hole natal kick, con- 
sidering five difTerent natal kick distributions (including a 
zero natal kick) and give this kick a random direction which 
we then add to the mass-loss kick due to the supernova to 
produce the total kick speed of the binary Vk- The gained 
velocity will be added randomly to the circular velocity of 
the binary within the Galaxy. Each system is then integrated 
forward in time within the Galactic potential for ~ 3 x 10^ 
yr (which is the MS-time of the 1.5 M© companion), and 
its position is noted at random times over the trajectory. 
We are thus able to produce an entire population of BH- 
LMXBs given the initial distribution of progenitor systems 
in the Galactic disc, their binary properties (separation and 
stellar masses) and the natal kick distributions for the black 
holes formed. 

We populate the disc of the Galaxy assuming the disc 
distribution of binaries to be proportional to the surface den- 
sity of stars S(i?) ~ Soe 
from the Galactic centre of Rmax 
be the length scale of the thin disc of the Galaxy, where the 
progenitor systems are thought to be produced, Rd ~ 2.6 
kpc (McMillan 2011). Concerning the z-distribution of the 
binaries, we model it as an exp onential with scale height 
~ 0.167 kpc (JBinnev et al.lll988l ). 

The population is formed by 100 binaries with the fol- 



^'^'', with a maximum distance 
10 kpc. We chose Rd to 



lowing parameters: Mi = 11 Mq, M2 = 1.5 Mq, Mbh = 
7.8 M0, a = 10 R© (Ml is the mass of the progenitor of 
the black hole, M2 is the companion star and a is the pre- 
SN orbital separation). For the black hole mass, we choose 
the average ma ss of stellar black holes in the Galaxy (see 
lOzel et al.ll201(J ): for the initial orbital separation, our choice 



Figure 5. Binary population synthesis for a sample of BH- 
LMXBs. No natal kick has been imparted to the BH. Smaller 
dots correspond to the synthetic population, bigger ones to the 
observed binaries and the position of the Sun is denoted with a 
square. 



is guided by the typical results of common envelope evolu- 
tion considered in detailed binary evolution calculations for 
progenitor systems. We choose a typic al mass l oss in the su- 
pernova explosion of ~ 3 Mq (see iFrver fc Kalogerai i 20011). 
which delivers an associated mass-loss kick of ~ 20 km/s. 

In addition, the black hole receives a natal kick drawn 
from one of the following five distributions: 1) a natal kick 
of zero km/s ; 2) one drawn from the Hansen & Phinney 
distribution (jHansen fc Phinnev 1997); 3) one drawn from 
the bimodal distribution of lArzoumanian et al.ll2002l : 4) as 
2) but with the kick speed multiplied by the factor Mna/Mbh; 
and 5) as 3) but with the kick speed multiplied by the factor 

Mns/Mbh. 

We plot the positions of the 100 binaries -at random 
times of the trajectory- in Galactic cylindrical coordinates 
for zero black- hole natal kicks in Figure [5] and for the other 
four natal kick distributions in Figure |S] From Figure [S] it 
is clear that it is impossible to place BH-LMXBs seen at 
larger values of z when the black holes receive zero natal 
kicks. Either the Hansen & Phinney or the Arzoumanian 
distributions appear to fit the observed distribution of BH- 
LMXBs, whereas those produced by natal kick distributions 
with velocities reduced by a factor of Mbh/Mns (bottom pan- 
els in Figure [Sjl appear to produce distributions which are 
more concentrated on lower values of z. 



5.1 Statistics of the results 

In order to quantify the results of the BPS, we show in 
figure [7] the fraction of binaries that at some time over the 
trajectory are located at a distance z from the Galactic plane 
less than a certain value. We include in the plot the results 
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Figure 6. Binary population synthesis for a sample of BH-LMXBs. Natal kicks have been drawn from Hansen & Phinney (top left) 
distribution, a bimodal distribution (top right), whereas the bottom figures correspond to the reduced natal kicks. Smaller dots correspond 
to the synthetic population, bigger ones to the observed binaries and the position of the Sun is denoted with a square. 



of the BPS for which no NK has been imparted to the black 
hole. 

It is evident how the mass-loss kicks alone cannot ac- 
count for the z-distribution of the observed binaries. A re- 
duced Hansen & Phinney NK cannot make the binaries that 
are located at z > 1 kpc; in particular, the percentage of 
binaries that get to z higher than 1 kpc is only the 0.5%. 
With a reduced Arzoumanian NK the percentage gets higher 
(~ 6%), though the fit remains unsatisfactory (see table |4]). 
It then turns out to be very difficult, with a momentum 
conserving kick, to reproduce the binaries XTEJ118-I-480, 
1705-250, XTEJ1859+226, which are located respectively at 
z = 1.5, 1.36, 1.20 kpc. Section[6]is dedicated to the detailed 
study of these sources. 

We are aware that the integration time we chose (~ 
3 X 10^ years) might be higher than the actual age of some of 
the observed binaries, particularly those whose mass trans- 



fer is driven by angular momentum losses, or those with 
a relatively massi ve companion star (M 2 ~ 3 Mq, see the 
updated catalogue iRitter fc Kolbll2003l ). We then carry out 
other two synthesises for the BH-LMXBs, integrating their 
trajectories for ~ 10* years and for ~ 5 x 10® years. In 
the first integration, the percentage of binaries that reach z 
higher than 1 kpc is 0% for a reduced Hansen & Phinney 
NK and 1% for an Arzoumanian reduced NK. In the sec- 
ond integration, the percentages are respectively are 0.2% 
and 2.2%. Concerning the KS-test, the resulting probabili- 
ties get 1 or 2 order of magnitudes lower when choosing a 
reduced integration time; this is easily explained, since the 
binary doesn't live long enough to be seen at high z. 

We wonder whether a larger mass-loss kick would af- 
fect our conclusions. Referring to equation [6l we see that 
the mass-loss kick increases either in the case of a larger 
mass-loss, or a more compact initial binary, or a larger com- 
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Figure 7. Cumulatives which show the fraction of BH-LMXBs 
versus the distance from the Galactic plane, for the four differ- 
ent natal kicks (dotted line is for an Hansen & Phinney NK, 
dotted-dashed line is for a reduced Hansen & Phinney NK, short- 
long-dashed is for a reduced Arzoumanian NK, whereas a zero 
NK scenario corresponds to the long-dashed line) . For the Arzou- 
manian NK, we tested two additional scenarios, a NK lying in the 
orbital plane and a NK perpendicular to it (see the three short- 
dashed lines). Cumulatives are to be compared with the observed 
one (solid line). 



Figure 8. Cumulatives which show the fraction of NS-LMXBs 
versus the distance from the Galactic plane, for the four different 
natal kicks (dotted line is for an Hansen &: Phinney natal kick, 
short-dashed line is for an Arzoumanian natal kick, whereas a 
zero natal kick scenario corresponds to the long-dashed line). Cu- 
mulatives are to be compared with the observed one (solid line). 



Table 4. KS probabilities for the EPS. 



panion mass. It is believed that in the SN event the Hehum 
star loses no more than 3 — 4 M© (before exploding as a 
SN, the Helium star suffers from strong WR winds after the 
comm on envelope phase, see for example iFrver fc Kalogeral 
I2OOII ). Concerning the initial orbital separation, there is a 
limiting minimum value for which either one or both of the 
two stars fill their Roche lobe. The following parameters. 
Ml = 11 Mq, M2 = 3.0 Mo, Mbh = 7.8 M©, a = 6 R©, 
give a recoil velocity Vmik of ~ 40 km/s. We then perform 
two synthesises in which we fix the mass loss kick to ~ 40 
km/s, testing the two types of reduced natal kicks. The cor- 
responding KS probabilities remain unsatisfactory: 2 x 10~^ 
for a reduced Hansen NK, and 2 x 10~^ for a reduced bi- 
modal NK. We shall also stress that the integration time 
for these two synthesises has been set to ~ 3 x 10® years; 
decreasing the integration time would make the KS proba- 
bilities even lower. 

We perform a BPS for NS-LMXBs as well (binary pa- 
rameters chosen: Mi = 3.5 Mq, M2 = 1.0 Mq, Mns = 
1.4 Mq, a = 7.0 Rq): in Figure [5] results are shown. It 
is pretty clear that a bimodal distribution better fits the 
observed sample. This a strong case for neutron stars re- 
ceiving a bimodal NK at birth. Previous studies, focused on 
N Ss in Be X-ray bina ries and double NS bi naries (see works 
bv lPfahl et "ai]|2002bl and lWong et al.ll2012l '). showed that at 
least some of the NSs should have received a lower kick at 
birth. We highlight our work as the first test of a bimodal 
distribution being a better fit to the Galactic position of NS 



KS probabilities 




BH-BPS 


NS-BPS 


Integration time 3 X 10® 


years 


Hansen NK 


0.20 


2.6 X 10-3 


Bimodal NK 


0.18 


0.78 


Hansen MCK 


2 X 10-3 


- 


Bimodal MCK 


1 X 10-2 


- 


zero NK 


2 X 10-* 


5 X 10® 


Integration time 5 X 10* 


years 



Hansen MCK 7 X 10-* 
Bimodal MCK 5 X 10-^ 

Integration time 10* years 
Hansen MCK 4 X lO"* 
Bimodal MCK 6 x lO"* 



low-mass X-ray binaries. When excluding from the test the 
observed NS binaries that are located at 2 > 2 kpc, the KS 
probability rises to 0.19; this is easily explained, since the 
change in normalization shifts the simulated curve towards 
the observed one. 

In table[4]KS probabilities for the different types of sce- 
nario are shown (the probabilities are reasonably a ccurate 
for our number of data points, see lPress et al.lll993l l. 
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6 DISCUSSION 

We now aim at deriving the minimum natal kick required 
to place the observed BH-LMXBs in their current locations. 
Rather than considering the properties of a general progen- 
itor binary system (i.e. stellar masses, black-hole mass, and 
binary separation), we use the observed properties for each 
known system, where possible, to more accurately calculate 
the mass-loss kick and the effect of any na tal kick on the 
particular system (see iRitter fc Kolbl 120031 for an updated 
catalogue of LMXBs in the Galaxy). 

For simplicity here, we consider a kick in the (optimal) 
direction perpendicular to the Galactic disc. We first com- 
pute, via conservation of energy, the minimum kick V± for 
the binary to reach the current position {Ro, z) after trav- 
eling in the Galactic potential, assuming the binary is born 
right over the Galactic plane at some radius Rq. Thus 



iVx2-f$(iio,0) = $(i?o,2) 



(8) 



where <& is the gravitational potential of the Galaxy. 

We show in Table [S] the minimum perpendicular kick 
V± and minimum natal kick Vnk required for those BH- 
LMXB systems which have relatively well-constrained bi- 
nary properties. One should understand that the natal kick 
value quoted in Table 5 is the absolute lower limit, assuming 
that the natal kick occurred in the perfectly optimal direc- 
tion. In practice, the natal kick required will be much larger 
in many realizations (i.e. different kick directions for the na- 
tal and mass-loss kicks). For example, for many directions, 
the necessary natal kick to reproduce the current location 
of XTE J 11 18-1-48 is up to 300 km/s (as was also seen in 
iFragos et ai]|2009l '). 

The minimum peculiar velocity is normally greater than 
80 km/s for binaries that are located ai z > 1 kpc, or for 
binaries that are closer in towards the Galactic center and 
at z > 0.6 kpc. Particularly, for systems that are located 
at -R < 3 kpc from the Galactic center, typical required 
velocities are greater than 100 km/s for the highest-z sys- 
tems. These high velocities cannot evidently be accounted 
for by a mass-loss kick alone. In other words, the cur- 
rent location of at least some systems clearly requires the 
presence of a black hole natal kick broadly in the range 
100 — 500 km/s. For other systems, the current location 
could be reached with the black hole having received no na- 
tal kick. We no te that our resul t s are c onsistent with those 
found earlier (INelemans et al.l 1 1999. IJonker &: Nelemani 



20091 ) 



2004IWillems et al.ll2005l . iDhawan et al.ll2007l . lFragos et ~ 



The system 1705-250 requires the largest minimum na- 
tal kick. This is because the system is located close to the 
Galactic centre and therefore has climbed out of a deeper 
potential well, assuming it was born at a comparable ra- 
dius in the Galactic disc. As there is a strong radial depen- 
dence of the Galactic potential this close in, we compute 
the minimum peculiar velocity both launching the binary 
out of the disc at _R = 0.5 kpc (the current location is at 
R = 0.55 kpc, z = 1.36 kpc) and R — 2 kpc. In the first 
case, we get V± ~ 400 km/s, while in the second case a ve- 
locity V± ~ 250 km/s is needed. These velocities require a 
minimum natal kick of 440 km/s and of 260 km/s respec- 
tively. Our results might have been affected by the choice of 
a spherically symmetric bulge (i.e. as = 0.0 kpc, see equa- 



tion [2] for the potential of the spheroid); we then take a 
pseudobulge with as — 1.0 kpc. The resulting V± is ~ 320 
km/s for R = 0.5 kpc and V± ~ 220 km/s for R = 2.0 kpc. 
The associated minimum natal kicks are 340 and 230 km/s 
respectively. For both of the two model for the Galactic po- 
tential, the required velocities are larger than the largest 
velocities drawn from a reduced-velocity kick, from either 
the Hansen & Phinney or Arzoumanian kick distributions. 
We may wonder whether our conclusions on the re- 
quired minimum natal kick would be affected by a new 
estimation of the distance. For all of the previously men- 
tioned four binaries, the distance has been derived from 
the estimation of the absolute magnitude of the compan- 



ion (see respectiv e ly Orosz et al. 2002 



Gelino et aij|2006l . iBarret et al.lll996l '). iJonker fc NelemansI 
20041 observed that this method typically underestimate 



Orosz et al.l 12001 



the distance. Also, we point out that any underestimation 
of the contribution of the disc to the observed magnitude 
of the companion star, would lead to an underestimation of 
the distance. To quantify the effect of a new estimation of 
BH binaries distance, we computed the minimum required 
natal kick, the distance being 10, 25, 50, 100 % larger than 
the nominal value. We perform the computation for our four 
candidates of BHs receiving the same NK as NSs. The re- 
quired minimum NK decreases in all the cases except for 
XTEJ1118+480. This is easily explained since a larger dis- 
tance would move the binary further out the Galactic poten- 
tial weU. The binaries 4U1543-47, 1819.3-2525, 1705-250 re- 
quire a minimum natal kick of 45, 22, 92 km/s, respectively, 
when the distance is increased of the 100%. XTEJ1118+480 
shows instead an increase of the minimum NK up to ^ 100 
km/s when the distance is multiplied by a factor of 2. 

An alternative scenario for the formation of BH bina- 
ries would be via dynamical interactions in globular clusters 
(GCs). However, it is still uncertain whether BHs are re- 
tained in GCs or whether they follow a different dynamical 
evolution than NS binaries. So far, no B H X-ray binary has 
been found in Galactic globular clusters (jVerbunt fc LewirJ 
120061 ): the stronge st BH candi date in a GC is the one 
found bvJMaccaro ne et al.ll2007l . in the Galaxy NGC 4472. 
The question whether black holes might be retained in 
a globular cl uster has been larg e ly discussed in the lit 
erature (e.g.. [ Kulkarni et al.l 119931, ISigurdsson fc Hernguist 



ll99AIPortegies Zw art fc McMillanil2000l . lMiller fc Hamiltoii 
2002j). It is generally thought that black holes would tend 
to decouple dynamically from the rest of the cluster and to 
segregate into the core, where they would form BH-BH bi- 
naries. Sequential dynamical interactions between these bi- 
naries and single BHs would lead to the ejection of the BHs 
from the cluster in a timescale shorter than ~ 10^ years. 
In case one black hole survives in the cluster, it could po- 
tentially capture a stellar companion via two main mecha- 
nisms: tidal capture of a star by the BH or e xchange interac- 
tions of the BH with a primordial binary (see'Kaloger a et al.l 
2004 and reference therein). Stars the BH is interacting with 
have a typical mass < 1 Mq at this stage of the life of 
the GC. It has been shown that the two-body tidal cap- 
ture scenario betwe en a NS and a low - mass star is likely to 
resul t in a merger ( JDavics et al."l992, Kumar fc GoodmanI 
Il996l . [McMillan etal...l990. . .Rasio fc Shapiro..l9S?II '). Specifi- 
cally, [Davies_et_alJ[l993 showed that an encounter between 
a NS and a red-giant would result in a merger in some 70% 
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of the cases, and some 50% for an encounter between a NS 
and a MS-star. For a BH, tidal forces are expected to be 
much larger, so that a merger becomes even more likely. Re- 
garding the exchange interaction scenario, we find it unlikely 
that the resulting BH-MS star binary will get a large kick 
in subsequent dynamical encounters for the binary to be ex- 
pelled from the cluster. Nevertheless, in the optimistic case 
that the binary managed to be ejected from the cluster with 
a velocity comparable to the escape speed of the cluster, we 
may compute the resulting overall distribution of Galactic 
BH binaries. We assume the binaries to be born in a spheroid 
of 20 kpc radi us around the Galactic c entre, taking the halo 
distribution of lDehnen fc Binnevlll99a . We then kick the bi- 
naries with a velocity of 45 km/s and we follow their motion 
in the Galactic potential for ~ 3 x 10^. The resulting KS test 
gives probabilities lower than 8 x 10~^, even when we double 
the BH binaries distance. This result is not surprising since 
we rarely get any binaries in the disc when assuming that all 
binaries are born in GCs. It could be that this mechanism 
would work for the BH binaries found at the highest z. How- 
ever, at least in the case of XTE J1118-I-480 there arc strong 
arguments for rejecting a GC origin. ^Gualandris et al. 2005 
estimated the age of the system, using stellar evolution cal- 
culations, to be betwee n 2 and 5 Gyr, renderi ng a globular 
cluster origin unlikely. [Hernandez et al.l 120081 performed a 
detailed chemical analysis of the optical star. Starting from 
different initial metallicities of the companion, they calcu- 
lated the expected abundances after contamination from SN 
nucleosynthesis products and were able to rule out a halo 
origin for this BH binary. Additionally, iFraeos et al.l l2009l 
claimed that the surface metallicity of the donor star right 
before the onset of RLO might have been even higher than 
the observed one, which would make the argument for a disc 
origin stronger. A natal kick seems to be required for this 
system and one which exceeds the range of kicks obtained 
from either the Hansen & Phinney or Arzoumanian distribu- 
tions with a reduction by the ratio of black hole to neutron 
star masses. 

From table [5] we see that natal kicks exceeding 70 km/s 
are required for several systems: 4U1543-47, 1819.3-2525, 
XTEJ1118-I-480 and 1705-250. These binaries provide us 
with evidence that black holes receive natal kicks of the same 
size as those received by neutron stars. One might have ex- 
pected the black hole natal kicks to carry the same momen- 
tum as those for neutron stars if, for example, a neutron star 
formed first (and having received a kick) and then a black 
hole formed later as a result of fall-back material within the 
supernova. In particular, the magnitude of the natal kick im- 
parted to the BH depends on the competition between two 
timescales: the fall-back timescale r/6 and the timescale of 
the mechanism leading to the natal kick r„fc . If r/f, > Tnk , we 
expect the fall-back material not to receive the same natal 
kick has the proto neutron star; the BH natal kick will then 
be reduced by the ratio Mns/A^bh- In case r/i, < T„k, we 
expect the BH to receive a full natal kick. From our simula- 
tions it seems that at least in some of the cases the fall-back 
material received the same NK as the proto-NS. Our result, 
already strongly suggested in the overall distribution of z 
seen in the binary population synthesis, is surprising. How- 
ever, one should note the large natal kick speeds obtained 
by the bla ck hole in some of the supernova simulations per- 
formed bv lFragos et~ai]|2009l . 



We would also like to point out the recent measure- 
ment of the distanc e of the BH candidate MAXI J1659-152 
(jJonker et al.l 120121 '). Taking the medium value for the dis- 
tance in the allowed range d = 6 ± 2 kpc, we derive a dis- 
tance from the Galactic plane of z ~ 1.7 kpc. This candidate 
would add to the number of BH binaries fou nd at large dis- 
tance from the plane of the Galaxy (see also iKuulkers et al.l 
[2012), thus likely enlarging the sample of strong candidates 
for BHs receiving large kicks. 

A number of mechanisms have been suggested for na- 
tal kicks received by neutron stars, some of which will 
also apply to black holes, especially those forming through 
the subsequent fall-back of material onto a proto neutron 
star. Suggested kick mechanisms include two main scenarios: 
hydrodynamically-driven kicks and neutrino-driven kicks. 
The former can either be caused by asymmetries in the con- 
vective motions under the stalled shoc k (sec Herant et al.l 
Il994 iBurrows et al.lll995l Ijanka fc Muelleij|l99^ . and the 
recent work bv lNordhaus et al.ll201ll ) or by over-stable osci l- 
lation modes of the progenitor core (|Goldreich et al .11 19961 1. 
The latter are produced by asymmetries in the neutrino flux 
in a strong magnetic field (Arras & Lai 1999, Lai & Oian| 
Il998h . Electromagnetically-driven kicks, instead, act once 
the neutron star has formed: t he off-centre rotating dipole 
impart the neu tron star a kick IjHarrison fc Tademarulll975l . 
[Laiet al. 2001). For a review of neutron star natal kicks, 
see lLai 2004. Alternatively, if the core is rotating extremely 
rapidly it may form a central object surrounded by a mas- 
sive disc on collapse, which may in turn fragment possi- 
bly producing a second compact object orbiting close to 
the central black hole or neutron star. This secondary will 
rapidly spiral-in towards the primary. If the secondary is 
a neutron star, then it may transfer mass to the primary 
until it reaches the minimum mass for a neutron star at 
whic h point it will explode potentially giving the p rimary a 
kick (|Colpi fc Wassermanll2002l . lDavies et al.]|2002bl ). In case 
both objects are black holes, then a merger kick may result 
from the asymmetric emission of gravit a tional waves (see 
iRosswog et al.1 |2000|. IZhuge et all Il994 iRasio fc Shapiro! 



1994, Ruffcrt fc Jankal 120011 ). In both cases, the kick re- 
ceived can be several hundreds km/s. Or else, when the 
disc on collapse forms directly around a BH, it might re- 
lease its energy as a powerful jet: if this jets happens to 
be one-sided, it might impart a kick to th e central BH, 
as suggested by [Barkov fc Komissaroy| l20ld . The BH-BH 
merger scenario might be tested via BH spin measurement. 
Evidence of highly rotating BHs comes from the properties 
of X-ray spectra of Galactic BH binar ies (see for exampl e 
IZhang et al]|l997l . lLaoilll99ll and also I Fender et al.ll2010l '). 
Ruling out the accretion of matte r from a low-mass c ompan- 
ion as origin of the BH spin (see lMcClintockll2006|). a BH- 
BH m erger might well fit in this scenario. I Herrmann et al.l 
120071 estimated the recoil velocity for two coalescing BHs of 
equal mass and opposite and equal spin. The recoil velocity 
reaches ~ 470 km/s in the extreme case. iMerritt et al.ll2004l 
considered the case of two BHs of unequal mass; they cal- 
culated a maximum recoil velocity of 450 km/s for a mass 
ratio q in the range 0.2 — 0.4. 
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Table 5. V± necessary to get to the observed position, and cor- 
responding minimum natal kick Vnk "■ 



Name 


Vx 


Vnk 


R 


z 




(km/s) 


(km/s) 


(kpc) 


(kpc) 


3 4U 1543-47 


95 


80 


3.92 


0.70 


XTEJ1550-564 


22 


10 


5.0 


-0.14 


GROJ1655-40 


36 





4.98 


0.13 


1659-487 


113 


- 


3.25 


-0.67 


1819.3-2525 


160 


190 


2.14 


-0.82 


GRS1915-I-105 


5 





6.62 


-0.03 


GS2023-I-338 


10 





7.65 


-0.09 


GROJ0422-(-32 


25 


10 


9.91 


-0.41 


A0620-003 


10 





8.92 


-0.12 


GRS1009-45 


40 


15 


8.48 


0.62 


XTEJ1118-I-480 


80 


70 


8.73 


1.50 


1124-683 


50 


40 


7.63 


-0.73 


XTEJ1650-500 


20 


- 


5.71 


-0.15 


1705-250 


420 


450 


0.55 


1.36 


XTEJ1859-I-226 


80 


- 


7.23 


1.20 


GS2000-I-251 


15 





7.21 


-0.14 



" For BH-LMXBs that lack strong observational constraints, we 
are unable to calculate accurately Vnk so leave it blank here. 



7 CONCLUSIONS 

In this Paper, we have considered the distribution of low- 
mass X-ray binaries containing black holes (BH-LMXBs) 
within the Galaxy as a function of the distribution of natal 
kicks given to the black holes. We have synthesized a BH- 
LMXB population by forming systems randomly through- 
out the Galactic citealt, weighted by stellar surface den- 
sity of the disc. We have given each binary a mass-loss 
kick due to the supernova explosion of the primary and 
added to this kick a black hole natal kick drawn from one 
of five kick distributions; 1) a natal kick of zero km/s; 

2) one drawn from the Hansen & Phinney distribution 
IjHansen fc Phinnev 1997); 3) one drawn from the bimodal 
distribution of I Arzoumanian et al.l 12002 : 4) as 2) but with 
the kick speed multiplied by the factor Mns/Mhh', and 5) as 

3) but with the kick speed multiplied by the factor Mns/Mbh- 
We have added the two kicks together with random direc- 
tions and combined (randomly) with the original orbital ve- 
locity within the Galaxy. The trajectory of each binary has 
then been integrated within the Galaxy. 

A number of observed BH-LMXBs are found in excess 
of 1 kpc from the Galactic disc. By comparing our synthe- 
sized population to the observed systems, we show that the 
hypothesis that black holes only rarely receive a natal kick 
is ruled out at very high significance. The computed dis- 
tribution is most similar to the observed distribution when 
the black hole natal kicks are drawn from the same velocity 
distribution as for neutron stars. Although we are unable to 
rule out that black holes receive smaller kicks than neutron 
stars, in a number of cases the required natal kick is very 
likely to exceed the maximum possible kicks in the reduced- 
velocity distributions (i.e. distributions 4 and 5 above). 
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